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Abstract 

We use the method of T-convergence to study the behavior of the Landau- 
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de Gennes model for a nematic liquid crystalline film in the limit of vanishing 
thickness. In this asymptotic regime, surface energy plays a greater role and 
we take particular care in understanding its influence on the structure of the 
minimizers of the derived two-dimensional energy. We assume general weak 
anchoring conditions on the top and the bottom surfaces of the film and the 
strong Dirichlet boundary conditions on the lateral boundary of the film. The 
constants in the weak anchoring conditions are chosen so as to enforce that a 
surface-energy-minimizing nematic Q-tensor has the normal to the film as one 
of its eigenvectors. We establish a general convergence result and then discuss 
the limiting problem in several parameter regimes. 


1 Introduction 

In this paper we study equilibrium configurations of thin nematic liquid crystalline 
films. Nematic materials are typically composed of rod- or disk-like molecules and 
can flow like fluids, yet they retain a degree of molecular orientational order similar to 
crystalline solids. There are several mathematical frameworks to study the nematics, 
leading to different, but related variational models. The focus of the present work 
is on rigorous dimensional reduction of the Landau-de Gennes Q-tensor model from 
three to two dimensions. 

We begin by briefly reviewing the basic theory of nematics. The local orientations 
of molecules in uniaxial nematics are described by a director—a unit vector in a 
direction preferred by the molecules at a given point. The director field underlies 
the Oseen-Frank theory [I] based on an energy penalizing for spatial variations of 
the director. This theory incorporates various elastic modes (splay, bend, twist) 
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and interactions with electromagnetic fields and has generally been very successful in 
predicting equilibrium nematic configurations. However the Oseen-Frank approach 
is limited in that it prohibits certain types of topological defects, e.g., disclinations, 
as the constraint that the director must have a unit length is too rigid. A possible 
remedy was proposed by Ericksen [2j who introduced a scalar parameter intended to 
describe the quality of local molecular orientational order. 

Despite the fact that Ericksen’s theory is capable of handling line defects, it still 
assumes that a preferred direction is specified by the director, excluding a possibility 
that the nematic can be biaxial. Here a biaxial state differs from a uniaxial state in 
that it has no rotational symmetry; instead it possesses reflection symmetries with 
respect to each of three orthogonal axes (only two of which need to be specified). 
Biaxial configurations are conjectured to exist, e.g., at the core of a nematic defect 
[3]. Further, certain nematic configurations cannot even be orientable, that is, they 
cannot be described by a continuous director field [4j. These deficiencies can be 
circumvented within the Landau-de Gennes theory in three dimensions that we will 
discuss in subsequent sections (see also H. 0 , and ®). Briefly, the Landau-de Gennes 
theory is based on the Q-tensor order parameter field that is related to the second 
moment of local orientational probability distribution. The relevant variational model 
involves minimization of an energy functional consisting of elastic, bulk and surface 
contributions. Recently there has been considerable activity on modeling with surface 
energy terms using both Q-tensor and director theories, including for example, 00 

0 0 unun]. 

We are interested in proper reduction of the Landau-de Gennes model to two 
dimensions in the thin film limit. In this asymptotic regime, surface energy plays a 
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greater role and we take particular care in understanding its influence on the structure 
of the minimizers of the derived two-dimensional energy. To achieve this goal we 
employ the tool of T-convergence that has proved successful in tackling problems of 
dimension reduction in other settings, such as elasticity [12] and Ginzburg-Landau 
theory [13]. We work in the domain x (0 ,h) where 0 < h <C 1 and 12 C M 2 is 
bounded and Lipschitz. In a subsequent publication we extend this analysis to the 
case of a small neighborhood of an arbitrary smooth surface, either with or without 
boundary, as a rigorous analog of the dimension reduction procedure in na (see also 

|I5|)- 

In Section [3] we introduce and analyze the general expression for the surface energy 
and then, in Section [I] combine it with the bulk and elastic terms to form the full 
non-dimensionalized three dimensional energy functional. In Section [5] we derive the 
expression for the limiting energy Fq and in Section [6] we analyze minimizers of Fo 
in different parameter regimes. In all regimes that we consider, it is crucial that the 
space of competing Q-tensors is constrained to accommodate the requirement that 
each tensor has a normal to the surface of the film as an eigenvector. This condition 
is forced in the limit by the surface energy and justifies the reduced Q-tensor ansatz 
imposed in [16], .17] in relation to experiments in [T8] . 

2 The Q-tensor 

In the three-dimensional setting, one describes the nematic liquid crystal by a 2-tensor 
Q which takes the form of a 3 x 3 symmetric, traceless matrix. Here Q(x ) models 
the second moment of the orientational distribution of the rod-like molecules near x. 
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The tensor Q has three real eigenvalues satisfying Ai + A 2 + A 3 = 0 and a mutually 
orthonormal eigenframe {1, m, n}. 

Suppose that Ai = A 2 = — A 3 / 2 . Then the liquid crystal is in a uniaxial nematic 
state and 


Q = ——1 <8) 1 — — m ( 
v 2 2 


m + A311 <8) n = S ( n 


n- -I 

3 


( 1 ) 


where S := ^ is the uniaxial nematic order parameter and n 6 § 2 is the nematic 
director. If there are no repeated eigenvalues, the liquid crystal is in a biaxial nematic 
state and 


Q — Ail (8) 1 + A 3 n <g) n — (Ai + A 3 ) (I — 1 <8> 1 — n <8> n) 

= Si (l ® 1 - 

where Si := 2Ai + A 3 and S '2 = \i + 2 A 3 are biaxial order parameters. Note that 
uniaxiality can also be described in terms of Si and S' 2 , that is one of the following 
three cases occurs: Si = 0 but S'2 7^ 0, S 2 = 0 but Si 7^ 0 or Si = S 2 7^ 0. When 
Si = S '2 = 0 so that Q = 0 the nematic liquid crystal is said to be in an isotropic 
state associated, for instance, with a high temperature regime. 

From the modeling perspective it turns out that the eigenvalues of Q must satisfy 
the constraints mm- 



A i G [-1/3, 2/3], for i = 1,2,3. 


(3) 
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3 Surface energy 


In this section we discuss the behavior of the nematic on the boundary of the sample. 
Here two alternatives are possible. First, the Dirichlet boundary conditions on Q 
are referred to as strong anchoring conditions in the physics literature: they impose 
specific preferred orientations on nematic molecules on surfaces bounding the liquid 
crystal. In the sequel we impose these conditions on the lateral part of the cylindrical 
sample dfl x (0, h ). An alternative is to specify the surface energy on the boundary 
of the sample; then orientations of the molecules on the boundary are determined as 
a part of the minimization procedure. We adopt this approach, referred to as weak 
anchoring , on the top and the bottom surfaces of the him. 

We seek a general expression for surface energy that has a family of surface- 
energy-minimizing tensors with the normal to the surface of the liquid crystal as 
their eigenvector. The requirement that the normal to the him is also an eigenvector 
of the Q tensor is motivated by the desire to model both homeotropic and parallel 
anchoring—corresponding to the nematic molecules oriented perpendicular and par¬ 
allel to the surface of the him, respectively p]. In the former case the uniaxial nematic 
tensor is prescribed on the boundary (up to the multiplicative order parameter) with 
the director being perpendicular to the surface of the him. In the latter case, the 
director orientation is perpendicular to the normal to the him but otherwise may be 
arbitrary. 

Consider the ’’bare” surface energy (Eq. 7 in [20] . see also [21] and [22] ) 

f s (Q, v) '■= Ci(Qu ■ u) + c 2 Q ■ Q + c 3 (Qv ■ u ) 2 + c a \Qv\ 2 , (4) 
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where q, i — 1,... ,4 are constants, A ■ B = tr (B t A) for any two n x n matrices A 
and B , and z/ G § 2 is a normal to the surface of the liquid crystal. This expression can 
be supplemented, in principle, by the surface Landau-de Gennes-type expression (Eq. 
75 in [ 6 ]) in order to control eigenvalues and to relax constraints on the constants c, 
in © that will be imposed below. This would amount essentially to augmenting (3D 
with an expression of the form (1T5D below. 

Fix v G S 1 and let 

A ■— {Q € : trQ — 0} . (5) 

We now explore different parameter regimes associated with minimization of the 
surface energy over the set A. Some comments along these lines can be found in [23| . 

Theorem 3.1. The minimum of f s (Q,v) over Q G A is achieved in the following 
five cases as characterized below in terms of the parameters c 4 ,..., c 4 : 

(i) //min{c 2 , 2 c 2 + c 4 , 3c 2 + 2c 3 + 2c 4 } > 0, then the minimum of f s (Q,u) is achieved 
at any Q that is uniaxial and homeotropic with the eigenvalue — , V>+0 C ^ cor¬ 
responding to the eigenvector u. 

(ii) If min{c 2 ,2c 2 + c 4 } > 0 and 3c 2 + 2c 3 + 2c 4 = c 4 = 0, then the minimum of 
fs(Q,v) is achieved at any Q that is uniaxial and homeotropic, i.e., with the 
nematic director parallel to v. 

(Hi) If min{c 2 ,2c 3 — c 2 } > 0 and 2c 2 + c 4 = 0, then the minimum of f s (Q,v) is 
achieved at any Q that satisfies Qu ■ v = an d has one eigenvector orthog¬ 

onal to v with eigenvalue 4c3 c f 2c o ■ 

(iv) //min{ 2 c 2 +c 4 , c 3 +c 4 } > 0 and c 2 = 0, then the minimum of f s (Q, u) is achieved 
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at any uniaxial or biaxial Q of the form 


<? = „ ( 6 ) 

where {m, n} is an arbitrary orthonormal frame in the plane tangent to the 
surface of the liquid crystal, and p G R. is arbitrary. 

(v) If c 3 > 0 and c 2 = c 4 = 0, then the minimum of f s (Q,v) is achieved and any 
minimizing Q must satisfy Qv ■ v — — 

In all other cases, iiifgg^ f s {Q, v) = — oo. 

Proof. Without loss of generality we may assume that v = (0, 0,1) := z. In order to 
hnd Q G A that minimizes f s , observe that 

Q ■ Q = 2\Qu\ — {Qv ■ v) + Q 2 • Q 21 

where Q 2 € M'sfm i s a nonzero square block of (I — v ® v) Q (I — v <g) v). Then the 
expression for f s can be written as 

f s (Q, v) = cQx ■ v) + C 2 Q 2 ■ Q 2 + (c 3 - C 2 ) (x ■ v) 2 + (2c 2 + c 4 ) \x\ 2 , (7) 

where x := Qv G M 3 . The traceless condition for Q can be reformulated in terms of 
Q 2 and x as 

tr Q 2 + x ■ v = 0 . ( 8 ) 

Thus, we are looking for the minimum of (J7J) among all Q 2 G and iGK 3 that 

satisfy the condition (jHJ). 




Suppose first that c 2 < 0. Then, setting x = 0 and observing that 


/ 


a 0 


\ 


Q 2 



satisfies the constraint (jHJ) for any a E M, we conclude that inf x Q 2 f s = — oo. We 
leave the verification of other parameter regimes which result in mi x Q 2 f s = — oo to 
the reader and instead concentrate on the five cases laid out in the statement of the 
theorem. 

To this end, let c 2 > 0. Minimizing f s with respect to Q 2 subject to the constraint 
(jHJ) is equivalent to minimizing Q 2 ■ Q 2 subject to (JHJ). We have 


(9) 


2 Q 2 + AI 2 — 0, 


where A is the Lagrange multiplier. Combining (jHJ) and (|9l) , we find that A = (x ■ v) 
and 



( 10 ) 


Substituting this expression back into (J7J), we have 


inf f s = inf f s , 

A R3 


where 



(ii) 
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or, equivalently, 


fs(x, v) := Ci(x ■ u) + ( ^c 2 + c 3 + c 4 ) (x ■ u) 2 + (2c 2 + c 4 ) |(I - u ® v)x\ 2 , (12) 


Observe that if min { 3 c 2 + 2 c 3 + 2 c 4 , 2 c 2 + C 4 } < 0, then inf x f s = — 00 . Suppose now 
that min {3c 2 + 2c 3 + 2c 4 , 2c 2 + c 4 } > 0 so that there is a critical point x of f s that 
satisfies 

df s 


dx 


= c x u + (2c 3 - c 2 ) (x ■ u)u + 2 (2 c 2 + c 4 ) x = 0. 


It follows that x is parallel to u and 


x = 


Cl 


3 c 2 + 2c 3 + 2c 4 


-v 


(13) 


is a minimum. Consequently, 


Qv = Xu, 


(14) 


where A = 


minimum 


Cl 


3c2+2c3+2c4 ' 


Combining (fT4D with the expression for Q 2 , we find the single 


Q = y C ® v - li), 


(15) 


of the surface energy corresponding to a fixed uniaxial nematic state with the or¬ 
der parameter S = ^ and the nematic director u. This is the case of so-called 
homeotropic (perpendicular) anchoring. This establishes (i). 

Proceeding with the proof of (ii), if min{c 2 , 2c 2 + c 4 } > 0 and 3c 2 + 2c 3 + 2c 4 = 
Ci = 0 then the expression (1T21) for f s reduces to its last term. Hence minimization 
simply requires x to be parallel to u meaning that u is an eigenvector of Q. Further, 
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Q is uniaxial since the corresponding eigenvalue is x ■ v and Q 2 is given by (TTOl) . 

Next, if min{c 2 ,2c 3 — c 2 } > 0 and 2c 2 + C 4 = 0, then (fill) reveals that f s is a 
quadratic function of x ■ v that is minimized when x ■ v = C2 ^ 2ci =: Combining 
this observation with (TTUT) we obtain that the minimizing Q must be of the form 


1 -a /2 0 q 13 ^ 

Q = 0 —cr /2 q 23 


y q 13 Q23 cr 


From this we readily see that (—< 723 , < 713 , 0) is an eigenvector of Q. This establishes 
(iii). Note that the minimizing set of f s consists of a family of biaxial tensors that 
contains a single uniaxial representative corresponding to a homeotropic boundary 
condition. Here the biaxial tensors do not have a normal to the surface of the him as 
an eigenvector. 

Now we pursue the regime c 2 = 0. In this case, from ([TJ) we see that f s is 
independent of <52- Thus the minimizing matrix Q 2 is arbitrary as long as it satisfies 
the trace constraint (jSJ). If, in addition, min{2c 2 + C4, c 3 + C4} > 0 the energy is still 
minimized by x given by (fl3l) with c 2 = 0, i.e., 



It also follows that if {m, n} is an arbitrary orthonormal frame in the plane tangent 
to the surface of the liquid crystal, then any 
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minimizes the surface energy. This verifies (iv). 


Finally in case (v) the expression (J7J) reduces to a quadratic function of x ■ v that 
is minimized at x-v — —ci/2c 3 . Observe that in this case the minimizing Q 2 is again 
arbitrary up to the trace constraint ((HD - □ 

Remark 3.1. Note that the eigenvalues determined in the cases (i), (Hi), and (iv) 
must respect the bounds (J3]) on eigenvalues of Q thereby imposing additional restric¬ 
tions on the parameters ci,..., c 4 . 

Having explored all possible parameter regimes we now focus on case (iv) where 
min{ 2 c 2 + C 4 , C 3 + C 4 } > 0 and C 2 = 0. Here the degeneracy of the set of tensors 
minimizing the surface energy f s provides sufficient freedom for nontrivial reduction 
to two-dimensional limits that we will carry out in the next section. An alternative 
approach would be to extend the surface energy by including quartic terms [23J and 
even surface derivative terms m- 

By rearranging the terms in (Ill), the surface energy has the form 

f s (Q, v) = Ci{Qv ■ u) + c 3 (Qu ■ u) 2 + c A \Qu\ 2 

2 

+ c 4 |(I — u <g) u) Qu| 2 
= a [{Qu ■ v) - (3] 2 + 7 |(I - v <g> u) Qv | 2 , (16) 

up to an additive constant. Here a = C 3 + C 4 > 0, (3 = — 2 ( C .(+ C4 ) ■ an d 7 = c 4 > 0. 
This form of the surface energy explicitly demonstrates that the minimizer has v as 
one eigenvector with corresponding eigenvalue equal to (3. 


— ( c 3 + C4) 


(Qu ■ u) 


Cl 


2(03 + C4) 
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4 Landau-de Gennes Energy Functional. 
Non-dimensionalization 

Suppose that the bulk elastic energy density of a nematic liquid crystal is given by 

Ze(VQ) := —|VQ| H- —Qij,jQik,k H-~ Qik,jQij,k H- ~^~QlkQij,kQij,l (17) 

and that the bulk Landau-de Gennes energy density is 

hdciQ) ■■= atr (Q 2 ) + jtr (Q 3 ) + (tr (Q 2 ))\ (18) 

cf. [ 6 ], Here the coefficient a is temperature-dependent and in particular is negative 
for sufficiently low temperatures, and c > 0. One readily checks that the form (fj~8l) of 
this potential implies that in fact fmc depends only on the eigenvalues of Q , and due 
to the trace-free condition, therefore depends only on two eigenvalues. Equivalently, 
one can view fmo as a function of the two degrees of orientation Si and S 2 appearing 
in (J2J). Furthermore, its form guarantees that the isotropic state Q = 0 (or equiva¬ 
lently Si = S 2 = 0 ) yields a global minimum at high temperatures while a uniaxial 
state of the form (JT]) where either Si — 0, S 2 = 0 or Si = S 2 gives the minimum 
when temperature (i.e. the parameter a) is reduced below a certain critical value, cf. 
M. In this paper we fix the temperature to be low enough so that the minimizers 
of fLdG are uniaxial. We also remark for future use that fLdc is bounded from below 
and can be made nonnegative by adding an appropriate constant. In light of this, we 
will henceforth assume a minimum value of zero for fLdc- 

Let O C R 2 be a bounded domain with a Lipschitz boundary and let h > 0 be 
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Figure 1: Geometry of the problem. 


given (Figured]). Assume that the energy functional is 

E[Q] := [ {/ e (VQ) + f LdG (Q)} dV + [ f s (Q, z) dA , (19) 

JQx(0,h) Jflx{0,h} 

where z is a unit vector normal to the surface of the him. Given uniaxial data 
g G A) we prescribe the lateral boundary condition of the form 

Q(x, y, z ) = g(x, y) for {x, y ) G dVt and z G (0, h). (20) 

The admissible class of tensor-valued functions is then 

C 9 h := {Q G H 1 (G X (0, h)~ A) : QUx{ 2 } = 9, Vz G (0, h)} , 

where A is the set of three-by-three symmetric traceless matrices defined in (]5l) . 
Throughout this work we assume that g is uniaxial and is taken so that C 9 h is nonempty. 

It has been shown, however, in mi that when L 4 ^ 0 minimizers of (fT9l) may 
fail to exist. On the other hand, L 4 = 0 precludes an appropriate reduction to the 
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general form of the Oseen-Frank energy for nematics | 6 |. Since we are interested in a 
characterization of minimizers, we will set L 4 = 0 . 


We nondimensionalize the problem by scaling the spatial coordinates 


x y _ z 

x = —, y = —, z — —, 
D y D h' 


where D := diam(D). Set £ = 7 ^ 2 , e = and introduce f e (Q, VQ) := jf e (Q, VQ). 


Dropping tildes, we obtain 


/e(vQ) := [g im , jQim,j -|- M2Qik,kQij,j “ 1 “ -^3 Qij,kQik,j 


- [M2Qij i jQis i s + MsQisjQij^] 


H - 2 zQim ,3 + (M 2 + M 3 )Q i 3 i 3Q i 3 )3 ], (21) 


where M 2 = M 3 = the indices i,m — 1,2,3, and j,k,l = 1,2. Rescaling the 
Landau-de Gennes potential via fj^ciQ) ^ fLdc{Q ) an d dropping tildes gives 

f LdG (Q) = 2A tv (Q 2 ) + tr (Q 3 ) + (tr (Q 2 )) 2 , ( 22 ) 


where A — “, B — and iuj = Letting a = f, 7 = 7 , setting 


fs(Q,z) ■= fi 2 Q£fs(Q,z), 


and dropping tildes, the expression for the nondimensionalized surface energy is 


f s (Q,z) = ol [(Qz ■ z) - / 3] 2 + 7 |(I - z®z) Qz\ 


(23) 
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We conclude that the total energy is 


E[Q\ = £D 2 h [ (/ e (VQ) + Wl f LdG m dV + £D 3 [ f s (Q, z) dA , (24) 

■Jnx( 0,1) isix{o,i} 

where the rescaled domain is denoted by the same letter Introducing the non- 
dimensional energy F e [Q\ := -AjE[Q], we hnd using (I24|) that 

FAQ] = [ (/e(VQ) + wJ LdG {Q )) dV + - [ f s (Q, z) dA. (25) 

i!Jx(0,l) e JSJx{0,l} 

5 Convergence of minimizers of F e [Q\ when e 0 

Assume that an appropriate constant has been added to the Landau-de Gennes energy 
to guarantee that F e [Q\ > 0 . We wish to consider a range of asymptotic regimes 
corresponding to different magnitudes of a and 7 . To this end, we will assume that 
a = «o + eot\ and 7 = 70 + £71 for some nonnegative constants cko, op, 7 o, 7 i- Then 
the surface energy density (f25l) can be written as 

f,(Q, i) = /f >(«, 2 ) + e/PW, I), (26) 

where 

/i° } : = «o i(Qz ■ z) - fd} 2 + 7o|(I - z ® z) Qz\ 2 , (27) 

and 

/i X) : = «i [{Qz ■ z) ~ A } 2 + 7i|(I - 5 (g) z) Qz \ 2 . ( 28 ) 
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As it will become evident below, we can assume that ckoM = 7 o7i = 0. Let 


fold := 


f„ |/,°(VQ) + w,f iia (Q) + 2 fP(Q, z)]dA if Q e H' g , 
+oo otherwise. 


(29) 


Here 


/e°(vg) 


Qim,jQim,j H - -^-2Qij,jQik,k H - -^3Qik,jQij,k^ 


the space 

H l g -.= {QeH\n-V)-.Q\ aa = g} 

and 

©:={Qe.4:/i 0 >(Q)=0}, 
for some boundary data g : dCl —> V. 

We now state our main theorem on dimension reduction via T-convergence. For 
those unfamiliar with the notion, we refer, for example, to [25]. Note that whenever 
necessary we will view H g as a subset of Cf via a trivial extension to three dimensions. 

Theorem 5.1. Fix g : dfl —>- T> such that F[ g is nonempty. Assume that —1 < M 3 < 
2, and — | — < M 2 . Then r-liny F e = F 0 weakly in C{. Furthermore, if a 

sequence {g £ } e>0 C Cf satisfies a uniform energy bound F e [Q e \ < C 0 then there is a 
subsequence weakly convergent in Cf to a map in F[ g . 

Proof. First, we demonstrate that one can always choose a trivial recovery sequence. 
Indeed, if Q e = Q G Cf\F[ g then liin e ^ 0 F e [Qf = +00 = F 0 [Q] and when Q e = Q e 
Hg then F e [Q e ) = F 0 [Q e ] = F 0 [Q\ for all e. 

For the lower semicontinuity part of the proof, consider an arbitrary sequence 
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{Qe} £> o C Cf converging weakly in Cf to some Q 0 G if*. It has been established 
in [26] (Lemma 4.2) and [24] that when the elastic constants satisfy the conditions 
— 1 < M 3 < 2 , and — | — < M 2 the integral of f e is weakly lower semicontinuous 

in if*(f 2 x ( 0 , 1 )) and, in fact, 


/ e (VQ) > C|VQ | 2 (30) 

pointwise for all admissible Q, where C > 0. Then using Sobolev embedding, one 
finds 

liminf F £ [Q £ \ > F 0 [Q 0 ]. 

£—>■0 

On the other hand, if {Q £ } £>0 C Cf converges weakly in Cf to some Q 0 G Cf/if*, 
then Q 0 depends on 0 and/or Q 0 is not V- valued. In the first case, invoking (l30lh we 
have 

liminf F £ \Q £ ] > dim inf — / IQezfdV = +00. 

e 2 7ox( 0 , 1 ) 

I 11 the second case, 

/ /j 0) (Qo,2 )<M*0 , 

Jflx{ 0,1} 

(cf. (TT 6 |) )>. thus by strong L 2 -convergence of traces on O x { 0 , 1 }, 

lirn - [ /i 0 ) (Q £ , z) dA = +00. 

e ^° £ ijlx{ 0 , 1 } 

Finally, since the uniform energy bound implies a uniform if*-bound, there exists a 
subsequence weakly convergent in if*(f 2 x ( 0 , 1 ); *4) to a limit Q 0 that is independent 
of z. Further, strong convergence of traces in L 2 implies that Q 0 G if*. □ 
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Remark 5.1. When M 2 = M 3 = 0, one can easily argue that the convergence of 
the subsequence is, in fact, strong. Indeed, F e [Q] = F 0 [Q] for every Q e Hg, hence 
F e [Q e \ < F 0 [Qo] for all e > 0. Since 


/ fL dG (Qe)dV+ / fW(Q e ,z)dA 

' Ox (0,1) Jttx{ 0,1} 


/ ( fLdciQo ) + 2/g 1 ^(Qo; z)) dA as e —» 0 


we have 


limsup f (\Qe,x\ F\Qe,y\ H- ^IQe.zl 2 ^) dV < f (|Qo,.t | 2 + \Qo,y\ 2 ) dA. 

e—>0 J Ox(0,l) V e / Jn 


Combining this with the lower semicontinuity of the L 2 -norm of the derivative, strong 
convergence in Cf follows. 


6 Minimizers of the T-limit in different regimes 

In this section we explore minimization of the T-limit for different parameter regimes. 
The added penalty terms that originate from the three-dimensional surface term have 
a potential to disconnect uniaxial states whenever the Landau-de Gennes term be¬ 
comes dominant. As we will show, this may result in formation of singular structures 
such as boundary layers of Allen-Cahn-type or vortices of Ginzburg-Landau-type with 
possible emergence of biaxiality (cf. [27]). 

In order to apply Theorem 15.11 we note first that a minimizer Q e of F f exists for 
every e > 0 by the direct method. Furthermore, selecting any G £ we observe 
that F e [Q e \ < F e [G\ = F 0 [G] and so we can extract a weakly convergent subsequence 
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{Qe}e> o such that Q £ — 1 Qo for some Q 0 in H l fj . Since by the basic properties of 
T-convergence any limit of minimizers is a minimizer of the limiting functional, we 
conclude that Qo minimizes F 0 . 

We will assume throughout this section that 12 is simply connected with a suffi¬ 
ciently smooth boundary. In what follows we will take 7 > 0 to be independent of e 
(i.e., 71 = 0 in (J2BD) and consider two distinct parameter regimes for a under several 
boundary conditions. 


6.1 Regime = 0 

First, let aq = 0 so that = 0 and /j 0) given by ([271) . Then V -valued maps have 
z as one eigenvector with corresponding eigenvalue (3. Note that these maps are not 
necessarily uniaxial, although they are required to be uniaxial on the boundary. 

There are two types of uniaxial D -valued maps: those corresponding to ji = —(3/2 
and /j, — (3 in ([ 6 ]), respectively. When /i = —(3/2 we have 


Q = 


3(3 

Y 



and the uniaxial Dirichlet condition is completely rigid as Q is equal to a constant. 
Alternatively, when ji — (3, one finds that 


Q 


-3P 



) 


where n is an arbitrary unit vector held on the plane and one has the freedom to 
choose uniaxial Dirichlet boundary data of any degree. Here by the degree of Q we 
understand the winding number of the planar S 1 -valued vector held e 2 ^ on <912 where 
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n = e*A Note that in this definition, it is assumed that Q G H l ^ 2 {dVl) but this does 
not preclude the possibility of the phase if} of n jumping by an odd multiple of 7r after 
one circulation around dfl. When this happens, the vector held n is discontinuous, 
but the held e 2 ^' is smooth. 

Case 1 . We begin by characterizing Q 0 in the hrst, topologically simpler case 

Q\dnx(p,i)=9’-=^p(z®z-^) ■ ( 31 ) 


Unless specihed otherwise, we find it preferable from this point on to use the following 
representation of Q G Hg invoked, for example, in nn and motivated by simulations 

in [T6] : 

^ Pi ~ § P 2 0 ^ 

Q = P 2 ~P\ — f 0 • (32) 

v 0 o p , 

It is a convenient change of variables in the setting when one eigenvector of the 
Q-tensor is parallel to the z-axis. Note that the boundary condition in the new 
coordinates is 

Plan = (pi,P 2 ) Ian = 0. (33) 


Applying the identity 


fe(Q) = l (2 + M 2 + M 3 ) |Vp| 2 + I (6 + M 2 + M 3 ) | V/5| 2 

Z O 

m 2 + m 3 

H- 2 -“ (PixPx - P^-yPy + PtyP* + frxPv) + \ M 2 + M a\ ( PixP2y ~ Pi y P2x ) (34) 
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from id, the expression for F 0 [Q] takes the form 

= h[p] := f {i|Vp| 2 + ^H-'(IpI) j dV. (35) 

where we have used that (3 is constant and that the integral of the Jacobian of 
p vanishes due to (l33]h Here p = (pi,p 2 ), the parameters M = 2 + M 2 + M 3 , 
S = \jM/wi , and 

W(t) = 41 4 + Ct 2 + D (36) 

with C = 6 (3 2 — 4B/3 + 4 A and D £ M. Note that C, in particular, varies with 
temperature through its dependence on the coefficient A. It is plausible to assume 
that C may change its sign in appropriate circumstances. 

One easily observes that the minimizer of (155|) . subject to the boundary condition 
(l33j) has a constant phase and so (j35j) reduces to a scalar minimization problem for 
the modulus p := |p|. The minimizers of (j35]l satisfy the Allen-Cahn type equation 

— A p + -! -W'(p) = 0 in n, p = 0 on dtt. (37) 

o 2 


The function p = 0 always solves this problem and, in fact, is the unique critical point 
and thus the minimizer if C > 0. Therefore the minimizing Q-tensor in this regime 
corresponds to a constant uniaxial state 


Qf3 


f 0 0 ^ 
0 -f 0 


\ o 0 (3 ] 


(38) 
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(cf. (132]) ) matching the boundary condition (j3T]i . When C < 0, however, the con¬ 
stant state Q = Qp looses stability once Ai(0) exceeds —2C/h 2 , where Ai(h2) is the 
first eigenvalue of the Laplacian. A minimizing nonconstant solution emerges in this 
parameter regime when the value of p enforced on dVt by the surface energy does not 
minimize the bulk Landau-de Gennes energy. Indeed we expect a boundary layer to 
form in the vicinity of dfl, bridging p = 0 to the minimum value of W in the bulk 
when 5 -C 1 see e.g. [28] . 

Further, the nonconstant minimizing configuration cannot be uniaxial everywhere 
in as can be seen from the measure of biaxiality introduced in [29] : 


where £(Q) 
have 


0 implies that Q is uniaxial. If we express £ in terms of p and (3, we 




/3 2 (V~/3 2 ) 2 

(V + 3/3 2 ) 3 ' 


(40) 


Since (3 is hxed, there are finitely many values of p where £(p, (3) vanishes. We conclude 
that if p is a nontrivial solution of the elliptic boundary value problem (1571) then the 
minimizer is necessarily biaxial almost everywhere. 


Case 2. Now we turn to the case of the boundary condition 


<5|aox(o,i) — 9 — 3/3 



(41) 


where n : dfl —> S 1 is such that n <g) n is smooth and may have a nonzero degree 
in the sense described above. Then the tensor Qo minimizes the energy F 0 given by 
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(J35]), subject to the boundary condition 


Q |an 


-3/3 



(42) 


Using the representation ()32lh we have that Qo can be represented by the vector p 0 
that minimizes F 0 and satishes 


P = -3/3 \ - ~, ni?r 2 J (43) 

on dfl where |p| = /f. In fact, F 0 is altered by an additive constant depending on g 
due to the presence of the null Lagrangian p\ x p 2 y — PiyPix hr (O that will not affect 
minimization. 

In order to better understand the behavior of Qo let 5 be small. Then, in general, 
when C in (foTH) is negative, the corresponding variational problem is of Ginzburg- 
Landau-type, but with a boundary layer bridging the equilibrium value yj—C/8 of 
p in the bulk to that of 3/3/2 enforced by the surface energy on <90. Furthermore, 
for topologically nontrivial boundary data for p, the minimizing vector held po has 
to vanish somewhere within a vortex core structure of a characteristic size 6 in f2. 
Recalling (1321) . we have that Po(^) = 0 corresponds to a uniaxial state at x with the 
director pointing along z , namely 

Qo(x) = ^ (z®z- ^ . 

In the case of C > 0—making W convex with a unique minimum at p — 0—we 
expect again that a boundary layer would form along dfl bridging the boundary value 
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of p, namely 3/3/2, to zero in the interior of 12. In all cases discussed in this section, 
we expect that symmetry breaking caused by surface energy will induce biaxiality 
within small sets where Qo experiences large variations. 

Remark 6.1. It is also possible to consider an intermediate asymptotic regime with 
non-zero degree Dirichlet data where minimality does not force biaxiality. Suppose the 
surface energy is taken to be substantially smaller than in those cases considered so 
far but where the relative strength of the Landau-de Gennes contribution is stronger. 
Suppose, for example, that 7 = £71 and a 0 = 07 = 0 for some 71 > 0 in expressions 
m , d 2 HD while say wi = ^ for some wo > 0 . In view of f[ 25 T) one sees that the new 
r -limit would now take the form 

_ J /„{/e°(Ve)+2 7 „|(I-z«l)Q7} dA if Q € Hy, 
ro[Q\ ■— < (44J 

+00 otherwise. 


where if/ := {Q G R 1 (12; A) : Q\aci = g , fiAciQ) = 0}. Here we assume the uniaxial 
data g takes the form 


g = s 


* 



for some planar vector field n where S* corresponds to the preferred value dictated by 
the requirement fidcid ) = 0 and suppose g has even degree so as to ensure that H l r] is 
nonempty. In this case, say for 7 0 large, one expects the uniaxial minimizer to simply 
undergo an out of plane rotation of its director in a neighborhood of the boundary, 
keeping the degree of orientation fixed, thereby smoothly bridging the boundary value 
n to z inside the domain so as to accommodate the cost of the remnant of the surface 
energy term residing in the T-limit. 
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6.2 Regime / S (1) (Q, z) = aq [(Qz ■ z) - (3} 2 

Now consider the case when a = ea\ and 7 = 70 for some 70 , ai > 0. Then 
fs°\Q, z) = 70 1(I — z <g) z) Qz\ 2 and = oq [(Qz ■ z) — (3] 2 and the set V 

consists of traceless symmetric tensors having z as one of its eigenvectors. When 
= 0 the asymptotic behavior of the limiting f un ctional Fq was considered in m 
as 8 —> 0. They characterize the mininrizer by identifying finitely many interacting 
vortices whose positions are determined via minimization of a renormalized energy 
of the type introduced in [30] for the Ginzburg-Landau functional. This description 
is very similar to the already discussed regime above which formally corresponds to 
= 00 , i.e., an eigenvalue corresponding to z is identically equal to f3 in hi. Hence 
we expect similar behavior for any finite value of a\. 

7 Structure of the singular region for topologically 
nontrivial boundary data 

In this section we provide further insight into the structure of the singular region 
that develops for topologically nontrivial boundary data in Case 2 of Section 16.11 
when C < 0 . 

Suppose for the moment that the dimensional reduction was carried out as in 
Theorem 15.11 for a sequence of functionals without the surface energy term. I 11 this 
case, the conclusion of the theorem would remain the same, with the exception that 
no additional constraints would have to be imposed on the target space of admissible 
tensors for the limiting problem. I 11 particular, any tensor field described by a field 
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of symmetric traceless 3x3 matrices in satisfying the appropriate boundary 

conditions on dfl would be admissible. 

Now recall that any tensor representing a uniaxial nematic as in (JU) can be asso¬ 
ciated with a point on MP 2 , i.e., on § 2 with antipodal points identified (since n and 
—n lead to the same Q-tensor). This means that the tensor held Q\da, corresponding 
to the Sprained boundary data (FI2D can be associated with one or more semi-circular 
arcs along the equator of § 2 . If the degree of p in (1431) is an even integer, then in terms 
of the image in MP 2 this corresponds to an even number of such half-equators. One 
can then smoothly deform the Q-tensor within 0 through a field of uniaxial nematics 
as indicated in Figure [2] to, say, a point in the interior of O where the director points 
either north or south. That is, the even number of semi-circles originally overlapping 
on the equator can gradually migrate towards the poles as a family of closed curves 
on MP 2 . 

Note that a sufficiently smooth uniaxial nematic tensor field can be described by 
a smooth director field [3] with the boundary data of degree equal to the half of 
that for Q\dn■ Hence the uniaxial nematic tensor held can be constructed even if the 
degree of the corresponding director field on the boundary is non-zero. Putting this 
another way, the director is § 2 -valued in Q and its § 1 -valued topologically nontrivial 
restriction to dfl can be bridged without forming a singularity in the interior of by 
the director “escaping into the third dimension”. 

If the degree of p in (1431) is an odd integer, then the image of in MP 2 corresponds 
to an odd number of half-equators of § 2 . This curve is not contractible in MP 2 , i.e., 
it cannot be smoothly deformed into a pair of points in § 2 . As Figure [3] illustrates, 
the shortest smooth deformation of an odd number of half-equators of § 2 is a single 
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Figure 2: A uniaxial nematic tensor Q with the director n can be identified with a 
projection matrix P = n 0 n, up to a translation and dilation. Let = Lfi(O) C R 2 
and let s be some parametrization of a circle of radius r G (0,1]. (a) The boundary 
data n|gn = e ls corresponds to an equator in § 2 ; (b) The boundary data for P 
corresponds to two half-equators in § 2 ; (c) For 0 < r 0 < 1, the half-circles migrate 
toward the north and south poles of § 2 , respectively; (d) At r = 0 half-circles contract 
to the poles. This corresponds to the director pointing up or down at r = 0. 

half-equator in § 2 , that is a closed geodesic in RR 2 . Then a uniaxial nematic Q-tensor 
field with non-contractible boundary data on dfl cannot be smooth everywhere in 
and must have a singularity at some x 6 Q as 0 will have to assume all values on a 
single half-equator in 8 2 in an arbitrarily small neighborhood of x. 

Staying with this broader class of admissible tensors for the moment, if we now 
let wi in (1291) be large but finite (this is equivalent to setting S in ([35]) small), any 
deviation of Q from a uniaxial nematic state-a minimum of fLd,G~ will result in large 
energy penalty. The smooth energy-minimizing tensor field will then be approxi¬ 
mately uniaxial nematic everywhere. That is, as described above, if the degree of 
Q\dn is even then there is no topological obstruction preventing Q 0 from being al¬ 
most uniaxial nematic everywhere in Since Qq has a fixed set of eigenvalues that 
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Figure 3: Behavior of Q for the boundary data of degree 3. The setup is the same as 
in the caption to Figure [2] and 0 < r\ -C 1. 

minimize fLdG and the director vector is an eigenvector of Q corresponding to the 
eigenvalue with the largest magnitude, the variations of Q throughout the domain 
can be interpreted as rotations of the eigenframe of Q. On the other hand, if the 
degree of be. the degree of p in (1431) . is an odd integer, then no such smooth 

uniaxial deformation exists going into O and there will be a small region of the size 
1 /■s/wl where the topological constraint will force Q 0 to be isotropic and/or biaxial 

(cf. 0EI). 

Let us now contrast this more relaxed target space with what occurs in our in¬ 
vestigation where the full energy functional contains the surface energy term. The 
effect of this is that the admissible tensor fields Q for the limiting problem are now 
P-valued. Since the normal to the surface of the him is then one of the eigenvec¬ 
tors of Q, there are only two types of possible uniaxial nematic states for admissible 
tensors—those with the director either perpendicular or parallel to the surface of the 
him (Figure S]). If 6 in (135]) is small, the potential term W would force the liquid 
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Figure 4: Geometry of the target space of uniaxial nematic T >-valued maps. Only the 
values along the equator and at the poles of § 2 are permitted. 

crystal to be uniaxial nematic throughout most of 0, thus making the director either 
perpendicular or parallel to the surface normal z. This orientational constraint is due 
to the strong influence of the surface energy and it makes director escape through 
the sequence of nematic states extremely costly. Given any topologically nontrivial 
planar boundary data for the director on dfl, we expect the director to remain planar 
as much as possible thus leading to formation of a localized region or regions with 
large gradients inside for any choice of topologically nontrivial pgn. 

At the core of any singular region that may develop, one would expect that the 
orientations of the molecules would be parallel to the surface of the him, while re¬ 
maining random within that surface. The orientational distribution function then 
is axially symmetric with respect to the normal to the him and the corresponding 
Q-tensor is uniaxial with the director being parallel to the normal (cf. |T5] ). Since 
only two types of uniaxial states are permitted for admissible tensors, the connection 
between the uniaxial states in the core and away from the core has to occur through 
a sequence of biaxial states via the mechanism of so-called eigenvalue exchange 
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8 Conclusions 


In this paper we have considered a rigorous dimension reduction procedure for nematic 
films described within the Landau-de Gennes model. We have established a general 
T-convergence result for the sequence of non-dimensional Landau-de Gennes energies 
as a small parameter £ related to the film thickness tends to zero. Although this 
result is applicable to any combination of coefficients in the expression ([4|) for the 
surface energy, we only consider cases when a Q-tensor minimizing f s (Q, v) in ([4j) has 
the normal to the surface of the film as one of its eigenvectors. 

In the thin-film limit, the dominant contributions to the energy come from the bulk 
elastic energy terms containing the derivatives in the direction normal to the surface 
of the him as well as from the surface energy terms that are independent of e. The 
energy penalty due to normal derivative terms forces the limiting energy-minimizing 
Q-tensor held to be independent of the spatial variable perpendicular to the surface 
of the him and reduces the domain of the problem from three to two dimensions. The 
presence of the surface energy terms that are independent of £ imposes constraints on 
the target space of admissible three-by-three Q-tensors. The resulting limiting bulk 
energy defined over the constrained set of Q-tensors consists of lower order bulk and 
surface energy contributions. 

Our results, along with those in HZ], demonstrate that the limiting problems 
for various parameter regimes can be studied by using the techniques developed for 
scalar Allen-Cahn type problems or Ginzburg-Landau-type vector-valued problems. 
Depending on the relationship between the coefficients of the surface and bulk energy 
terms, minimizing Q-tensor helds can develop boundary layers. For topologically 
nontrivial boundary data and large Wi in (l29lh the Q-tensor helds are characterized 
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by vortices with characteristic core size ~ 1 /y/w[. Note that even though 1 /y/w[ in this 
case is small, it is still much larger than £ and two-dimensional vortices correspond 
to disclination lines perpendicular to the surface of the him. 
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